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Nomenclature 


=  lift  force  derivative 
CL  =  Ci^qocS/mu 

Cip  =  roll  damping  coefficient 

Q,  =  C[pqxSd7/2u 

C„t  =  pitch  damping  derivative 

c;,  =  -  c^sd'/iiu 

=  damping  moment  derivative 
C;.  =  Cmiq<KSd2/2I  u 

Cnpm  =  Magnus  moment  derivative 

C'np,  =  C^gooSd>/2Iu 

d  =  aerodynamic  reference  diameter 

F  as  vector  force 

g  =  acceleration  due  to  gravity 

h  -  angular  moment  vector 

/  =  pitch  or  yaw  moment  of  inertia 

Ix  =  roll  moment  of  inertia 

k  =  pitch  or  yaw  radius  of  gyration 

L  =  lift  force 


L0  =  lift  force  derivative 

l  =  jet-damping  moment  arm;  I/mxs 

M  =  vehicle  mass 

M  -  vector  moment 

M  =  vehicle  mass  rate  of  change;  propellant  mass  flow  rate 

m  =  damper  mass;  missile  mass 

p  -  roll  rate 

Paii  ~  critical  roll  rate 

Pr  =  HP/2 

q  -  pitch  rate 

qoo  =  dynamic  pressure 

r  =  yaw  rate;  slag  mass  radial  position;  missile  c.g.  coning 

radius 

S  =  aerodynamic  reference  area 

T  =  rocket  thrust 

t  -  time 

u  =  missile  forward  velocity 

V  =  velocity  vector 


=  missile  lateral  velocity  due  to  coning 

x,  r,z 

=  nonrotating  coordinates 

X,  y,z 

=  body-fixed  coordinates 

Xs 

=  static  margin 

y 

=  center-of-mass  offset 

a 

=  angle  of  attack 

f 

=  angle  of  sideslip 

=  damping  ratio 

=  total  angle  of  attack 
=  zero  roll  rate  trim  angle  of  attack 
=  coning  angle  (Euler  angle) 
as  coning  angle  rate  of  change 
=  precession  damping  exponent;  p/paii 
=  precession  damping  exponent 

=  un 

=  damping  parameter 
=  yaw  moment  parameter 
=  processing  coordinates 
=  complex  angle  of  attack,  +  id 
as  trim  force 

=  spin  angle  relative  to  6  plane  (Euler  angle) 
=  spin  rate  relative  to  9  plane 
=  precession  angle  (Euler  angle) 

=  precession  rate 
=  precession  rates 
=  precession  rates 

=  (Pr2+^),/3 
=  natural  pitch  frequency 
=  angular  velocity  vector 


Introduction 

ANYONE  who  has  worked  to  some  degree  with  the  gyroscopic 
motion  of  spinning  missiles  or  spacecraft  has  undoubtedly  dis¬ 
covered  that  the  Achilles  heel  of  the  vehicle  with  regard  to  angle 
of  attack  or  attitude  stability  is  the  “yaw’*  moment  or  “side"  mo¬ 
ment  or  “Magnus-type”  moment,  all  of  which  refer  to  a  similar 
type  of  disturbance.  Because  the  general  epicyclic  angle-of-attack 
motion  of  a  spinning  missile  is  a  superposition  of  two  coning  mo¬ 
tions,  one  could  argue,  with  few  exceptions,  that  all  missile  dynamic 
instabilities  are  coning  instabilities  caused  by  a  yaw  moment  distur¬ 
bance  that  effectively  undamps  one  coning  motion  while  damping 
the  other.1"4  Nicolaides,  in  an  early  paper  on  the  flight  dynamics 
of  modem  ballistic  missiles,5  commented  that  “all  flight  dynamic 
instabilities  observed  by  the  writer  have  been  observed  to  be  ulti-  • 
mately  of  the  pure  circular  type  ”  Referring  to  “catastrophic  yaw” 
in  the  same  paper,  he  remarked  that  the  pitching  and  yawing  motion 
“begins  to  grow  and  may  very  soon  reach  extreme  values  such  that 
the  missile  in  flight  looli  more  like  a  propeller  than  an  arrow” 
Similarly,  spin-stabilized  spacecraft  have  exhibited  coning  insta¬ 
bilities  in  which  the  angle  between  the  spin  axis  and  the  space  fixed 
angular  momentum  vector  grows  to  excessive  values.  It  is  shown 
that  this  occurs  from  energy  dissipation  when  the  spin  axis  is  one  of 
minimum  moment  of  inertia.6  This  coning  motion  can  be  damped 
by  an  actively  driven  unbalanced  rotor  or  pendulum  mass  that  is 
quite  small  relative  to  the  mass  of  the  spacecraft7  A  similar  passive 
damper  consists  of  a  small  mass  of  viscous  fluid  contained  in  an  an¬ 
nulus  that  rotates  at  the  precession  frequency  and  damps  the  coning 
angle  when  the  spin  axis  is  that  of  maximum  moment  of  inertia.8 

In  both  cases,  the  action  of  the  small  mass  imbalance  produces  a 
small,  stabilizing  yaw  moment.  The  passively  damped  pendulum  or 
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annular  fluid  damper  produces  a  positive,  destabilizing  yaw  moment 
if  the  spin  axis  is  that  of  minimum  moment  of  inertia.  A  destabi¬ 
lizing  yaw  moment  is  also  produced  by  an  oscillating  liquid  in  a 
spinning  shell.  This  was  first  demonstrated  by  Stewartson  for  an  in- 
viscid  liquid  in  a  spinning  top.9  Yaw  moment  effects  similar  to  fluid 
sloshing  can  also  be  produced  by  internal  moving  parts  in  spinning 
projectiles.10-12 

Recent  launches  of  a  class  of  spin  stabilized  vehicles  consisting 
of  a  spacecraft  and  a  payload  assist  module  (PAM)  upper  stage  have 
exhibited  a  coning  instability  during  the  latter  part  of  the  rocket 
motor  bum.13-17  This  has  attracted  much  attention  in  the  engineer¬ 
ing  community  because  virtually  all  of  the  launches  consisting  of  a 
payload  attached  to  a  PAM  solid  rocket  motor  have  exhibited  this 
coning  growth,  whereas  no  other  comparable  vehicle  (with  perhaps 
one  exception13)  has  shown  such  behavior.  Various  attempts  have 
been  made  to  explain  the  source  of  the  observed  coning  in  terms  of 
a  compliant  nozzle,  nozzle  entrance  flow  asymmetry,  spinning  body 
flowfield  effects,  and  slag  sloshing  effects.  There  is  no  satisfactory 
explanation  to  date  of  the  actual  cause  of  the  coning  instability,  al¬ 
though  a  study  of  the  oscillatory  motion  of  a  small  spring-restrained 
mass  to  simulate  slag  sloshing  in  the  presence  of  thrust15  showed 
that  an  instability  could  occur.  The  analysis  led  to  the  not  surprising 
conclusion  that  the  oscillation  frequency  and  phase  must  be  such  to 
produce  a  yaw  moment. 

Because  of  the  striking  similarity  of  the  forcing  function  required 
to  cause  coning  instabilities  in  spinning  bodies,  whether  they  are 
spin-stabilized  spacecraft  in  a  near-moment-free  environment  or 
high-speed  missiles  subject  to  large  aerodynamic  forces  and  mo¬ 
ments,  it  is  instructive  to  derive  the  general  conditions  required 
for  spinning-body  coning  instabilities  and  apply  them  to  familiar 
missile  and  spacecraft  coning  phenomena.  It  is  convenient  for  this 
analysis  to  use  the  classical  Euler  angle  coordinates  because  of  the 
simplicity  afforded  by  the  conditions  of  quasisteady  coning  motion 
and  the  simple  geometrical  interpretation  of  the  disturbance  forces 
and  moments  required. 

Equations  of  Motion 

The  rigid-body  motion  of  a  nominally  axisymmetric  missile  or 
spacecraft  is  described  by  the  general  force  and  moment  equations 

-=[V]  +  wxV  (1) 

m 

A f  =  [i]  +  «xA  (2) 

where  us  is  the  angular  velocity  of  a  rotating  coordinate  frame  in 
space,  V  is  the  velocity  of  the  center  of  mass  relative  to  the  rotating 
coordinates,  h  is  the  angular  momentum  of  the  body  with  respect  to 
the  rotating  coordinates,  and  F  and  M  are  the  forces  and  moments 
acting  on  the  body  along  and  about  the  coordinate  axes,  respectively. 
The  classical  Euler  angle  coordinates  arc  shown  in  Fig.  1,  in  which 
X ,  F,  and  2  are  nonrotating,  space-fixed  axes,  and  £,  tj,  and  f  are 
the  rotating  coordinate  axes  from  which  we  derive  the  equations 
of  motion.  We  consider  the  axis  of  symmetty  of  our  nominally 
axisymmetric  body  to  lie  initially  along  the  space-fixed  X  axis.  The 
axes  £,  ij,  and  £  are  obtained  from  X,  K,  and  2  by  a  rotation  (or 
precession)  0  about  X  followed  by  a  rotation  6  (or  pitch)  about  r). 
The  angular  velocity  a?  of  our  rotating  frame  in  space  with 

respect  to  the  £,  77,  f  axes  is 

id  —  (0cos0,  0,  if  sin0)  (3) 

or,  for  small  0, 

id  %  (ir,  0,  00)  (4) 

If  we  consider  the  further  rotation  <p  about  the  coordinate  axis 
then  the  axes  x,  y,  z  are  conventional  body-fixed  axes  commonly 
used  to  describe  aircraft  or  missile  motion.  Note  that  the  body-fixed 
axes  y  and  z  rotate  about  f  with  the  spin  rate  0  with  respect  to  the 
££  plane  such  that  the  resultant  body  angular  rate  p  in  space  about 
£  or  x  is  the  coordinate  angular  rate  0  cos  0  plus  the  spin  rate  0,  or, 
for  small  0, 

P*0  +  0  (5) 


X 


Fig.  1  Classical  Euler  angle  coordinate  system. 


The  corresponding  pitch  and  yaw  rates  about  the  body-fixed  axes  y 
and  z,  from  Fig.  1,  are 

q  —  0cos0  +  00sin0  (6) 

r  =  00  cos0  —  0  sin0  (7) 

The  angular  momentum  vector  h  with  respect  to  the  rotating  coor¬ 
dinates  £,  q,  {  is,  for  small  0, 

h  —  (/j  p,  16,1 00)  (8) 

where  7X,  /,  and  I  are  mass  moments  of  inertia  about  the  axes  £,  t)t 
C  (orx,  y,  z),  respectively,  for  the  axisymmetric  body. 

The  moment  equations  of  motion  for  a  constant  mass  body  in  the 
Euler  angle  coordinates  are  obtained  by  substituting  Eqs.  (4)  and 
(8)  into  Eq.  (2),  which  gives 


-y  -  $  +  (VP  -  00)00  (10) 

Y  -  “(00)  +  (0~  W)0  (11) 

where  p  s  Ix/I. 

Yaw  Moment  Damping  or  Undamping 

Consider,  first,  the  moment-free  body  in  steady  coning  motion, 
0  =  const.  It  follows  from  Eqs.  (9-1 1)  and  Eq.  (5)  that 

p  =  const  (12) 

0  =  MP  03) 

0  =  p(l-/x)  (14) 

The  body  processes  about  X  at  the  rate  up  with  the  coning  angle  0 
and  spins  with  respect  to  the  ££  plane  with  the  rate  p{  1  —  p,).  Note 
that  for  p,  <  1  the  spin  rate  0  is  in  the  same  direction  as  the  coning 
rate  0  and  for  p  >  1  it  is  in  the  opposite  direction.  Yaw  moment 
undamping  is  the  growth  in  the  coning  angle  0  caused  by  a  positive 
net  yaw  moment  M;  in  Eq.  (1 1).  A  negative  yaw  moment  will  cause 
a  decay  in  0.  The  response  of  the  coning  angle  to  an  applied  yaw 
moment  Aff,  compared  with  the  response  to  a  pitch  moment  Af7, 
is  seen  most  readily  if  we  apply  both  pitch  and  yaw  moments  of 
the  form 


(15) 

11 

* 

(16) 

where  a  is  a  multiple  that  determines  the  size  of  the  pitch  moment 
relative  to  the  yaw  moment  If  we  substitute  Eqs.  (15)  and  (16)  into 
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Eqs.  (10)  and  (11),  we  find  for  the  quasisteady  conditions  9^0 
and  \j/  %  const  that  the  motion  is  composed  of  two  coning  motions 
with  coning  rates 


V'u  =  W~ 


OK 

HP' 


OK 

PP 


and  amplitudes,  from  Eq.  (1 1),  described  by 


6  -  kU20  =  0 


(17) 


(18) 


where 


A|,2  =  "  - -  —  ±k  f{ip  -  (19) 

2* -up  V  VP/ 

On  integrating  Eq.  ( 1 8),  we  see  that  the  two  motions  grow  and  decay 
according  to 


6  =  00  exp  (kL2t)  (20) 

For  a  pure  yaw  moment,  a  =  0,  the  coning  rate  is  unchanged  from 
Eq.  (13),  and  the  coning  angle  grows  exponentially  according  to 

0=e0cxp^-tSj  (21) 

Consider  a  spin-stabilized  spacecraft  with  pp  =  7T  rad/s  and  let 
k  =  0.02.  For  combined  pitch  and  yaw  moments  with  a  —  10  such 
that  the  pitch  moment  is  10  times  the  magnitude  of  the  yaw  moment, 
the  effect  of  the  pitch  moment  is  to  reduce  the  fast  coning  frequency, 
from  Eq.  (17),  by  only  about  2%  and  to  increase  the  magnitude  of 
the  exponent  A.ii2  by  about  4%.  The  fast  mode  =  0.98  tt  rad/s 
grows  exponentially  with  a  time  constant  r  =  1  /A.j  of  15 1  s,  and  the 
slow  mode  ^2  =  02/n  rad/s  decays  with  the  same  time  constant 

Spacecraft  Precession  Damper 
A  damper  that  has  been  used  to  damp  the  coning  motion  of  a  spin- 
stabilized  spacecraft  is  readily  described  in  terms  of  yaw  moment 
damping.  Consider  a  processing  spacecraft  that  contains  a  small 
pendulum  mass  free  to  rotate  about  the  spacecraft  axis  of  symmetry 
in  a  plane  a  distance  h  from  the  spacecraft  center  of  mass,  as  shown 
in  Fig.  2.  The  mass  could  be  a  slug  of  liquid  contained  in  an  annular 
cavity,  which  is  a  type  of  device  actually  flown  to  damp  spacecraft 
coning  motion.18"23  If  the  mass  or  slug  of  fluid  is  free  of  any  drag 
or  viscous  damping  forces,  it  will  assume  a  steady-state  orientation 
at  4>d  —  0  in  the  negative  £  direction  due  to  centrifugal  force  from 
the  steady  precession  if.  Because  the  spacecraft  spins  at  the  rate  4> 
relative  to  the  or  9  plane,  the  mass  rotates  at  the  rate  4>  relative 
to  the  spacecraft  Any  viscous  damping  or  drag  on  the  mass  will 
displace  it  at  some  positive  or  negative  angle  <f>d  from  the  9  plane, 
depending  on  the  direction  of  i> ,  giving  rise  to  a  yaw  moment  The 
yaw  moment  is  simply  the  component  of  centriftigal  force  rmjf2R 
acting  in  the  r)  direction  times  its  moment  arm  h ,  or 

Mk  =s  mi/r2Rh  sin  <J>  (22) 

where  the  magnitude  of  R  is  a  function  of  the  coning  angle  9.  For 
p  <  1,  the  spin  rate  </>,  from  Eq.  (14),  is  positive,  which  would 
produce  positive  as  shown  in  Fig.  2  and  a  positive  yaw  moment. 
In  such  a  vehicle,  a  passive  damped  rotor  or  viscous  fluid  annular 
damper  would  undamp  the  coning  motion.  Note  also  that  any  viscous 
or  drag  force  that  rotates  the  pendulum  out  of  the  $  plane  produces  an 
opposite  roll  torque  that  will  change  the  roll  rate  accordingly.  For  the 
slender  spacecraft,  /,  <  /,  yaw  moment  undamping  of  the  coning 
motion  is  accompanied  by  a  decrease  in  the  roll  rate.  For  spacecraft 
configurations  with  p  >  l ,  for  which  the  passive  precession  damper 
was  devised,  the  spin  rate  0  is  negative  for  positive  if,  and  the  mass 
will  be  driven  in  the  opposite  direction  such  to  produce  a  negative  <J> 
and  negative  yaw  moment.  Because  the  pendulum  rotates  relative  to 
the  spacecraft  in  the  same  direction  as  ir ,  in  this  case,  it  will  exert  a 
positive  roll  torque  on  the  spacecraft  via  the  viscous  or  drag  forces 
and  increase  its  roll  rate.  Precession  damping  can  be  accomplished 
in  the  slender  spacecraft,  p  <  1,  by  actively  driving  an  unbalanced 


rotor  with  a  torque  motor,  for  example,  to  produce  the  equivalent  of 
a  negative  <j>d  or  <J>  in  Fig.  2  with  a  positive  4>.  Such  devices  have 
been  developed.7 

Vehicle  Coning  with  Axial  Thrust  and  Variable  Mass 

The  variable  mass  accompanying  thrust  from  a  spin-stabilized 
rocket  motor  or  PAM  produces  a  destabilizing  effect  on  coning  due 
to  a  decrease  in  the  lateral  moment  of  inertia  and  a  stabilizing  effect 
from  jet  damping.  Some  authors  include  the  variation  of  moment  of 
inertia  in  their  definition  of  jet  damping.  The  net  effect  of  the  variable 
mass  appears  to  be  stabilizing  for  the  solid  rocket  motor/payload 
configurations  that  have  been  investigated.  The  rocket  thrust  is  a 
potential  source  of  destabilizing  yaw  moment  in  the  presence  of 
a  small  mass  imbalance  relative  to  the  processing  0  plane.  These 
effects  will  be  illustrated  with  somewhat  simplistic  models. 

Consider  the  solid  rocket  motor  shown  in  Fig.  3  in  steady  coning 
motion  and  assume,  for  simplicity,  that  the  propellant  bums  in  such 
a  manner  that  the  center  of  mass  remains  fixed.  The  total  lateral  rate 
is  if  9,  for  small  9,  and  jet  damping  of  this  lateral  rate  produces  an 
opposing  yaw  moment2* 


=  Mt2if9  (23) 

where  M  is  the  negative  rare  of  change  of  vehicle  mass  from 
expended  propellant  and  —M  is  the  propellant  mass  flow  rate. 
Similarly,  it  has  been  shown  that  a  jet  damping  moment  that  op¬ 
poses  the  rocket  spin,  due  to  spin  momentum  imparted  to  the  ejected 
propellant  gas,  is  of  the  form25 

Mk  =  MR2Fp  (24) 

where  R  is  some  maximum  radius  of  the  mass  of  ejected  gas,  as¬ 
sumed  to  be  in  the  form  of  a  disk,  and  F  is  a  function  of  the  radial 
variation  of  the  rotational  velocity  imparted  to  the  gas.  If  we  in¬ 
clude  variable  moments  of  inertia  in  the  angular  momentum  vector, 
Eq.  (8),  and  substitute  in  the  moment  equations,  Eq.  (2),  with  the  jet 
damping  moments,  Eqs.  (23)  and  (24),  we  obtain 


'  mr2f  L\ 
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where  we  have  assumed  for  quasi-steady  coning  motion  that  9  % 
S  **  0  and  if  const  in  the  pitch  equation.  Because  4//*,  ///, 
and  M  are  negative,  the  variable  mass  effect  on  moment  of  inertia 
is  destabilizing,  whereas  jet  damping  is  stabilizing.  If  we  make  the 
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further  assumption  that  the  radius  of  gyration  k  in  yaw  remains 
constant,  where  /  =  Mk2  and  M  is  the  vehicle  mass  at  any  instant, 
then  the  yaw  equation,  Eq.  (27),  can  be  written 

'-(Ms*- 0  08 

and  integrated  to  give 

(29) 


as  a  byproduct  of  the  propellant  combustion.  The  mass  would  have 
to  rotate  or  oscillate  relative  to  the  motor  case  such  to  maintain  a 
position  out  of  the  ££  or  9  plane  described  by  0  <  fa  <  180  deg 
in  Fig.  2  so  that  the  thrust  acceleration  along  £  acting  on  the  mass 
produces  a  positive  yaw  moment  M f .  If  the  motion  were  such  that 
the  mass  occupied  the  position  180  <  fa  <  360  deg,  a  negative  yaw 
moment  would  be  produced  that  would  damp  the  coning  motion. 
Mingori  and  Yam15  considered  a  mass  restrained  by  a  radial  spring 
oscillating  in  a  plane  perpendicular  to  £  and  located  a  distance  h  fore 
or  aft  of  the  center  of  mass,  similar  to  that  shown  in  Fig,  2  but  with  a 
spring  restraining  the  mass.  They  came  to  the  not  surprising  conclu¬ 
sion  that  the  mass  must  lie  out  of  the  plane  containing  the  transverse 
angular  velocity  vector  (the  ££  plane  of  Fig.  3  in  which  rj/9  is  the 
transverse  angular  velocity  about  the  £  axis),  such  to  produce  a  pos¬ 
itive  yaw  moment  There  is  no  satisfactory  physical  explanation  for 
the  spring-mass  system  parameters  and  the  frequency  required  for 
instability,  other  than  a  loose  analogy  with  fluid  sloshing.  More  re¬ 
cent  analyses  with  spherical  pendulum  models  for  the  slag  mass  also 
required  physically  unrealizable  parameters  to  match  the  observed 
PAM  motions.1617  By  comparison  with  the  passive  rotor  or  annular 
damper  of  Fig.  2,  a  small  mass  in  the  aft  end  of  the  motor  subjected 
to  centrifugal  force  from  precession  and  a  small  viscous  drag  force 
would  be  driven  out  of  the  ££  plane  such  to  cause  a  negative  yaw 
moment  due  to  thrust  that  would  damp  the  precession  for  a  prolate 
vehicle  (/*  <  1).  Nevertheless,  the  amount  of  mass  imbalance  re¬ 
quired  to  cause  instability,  if  it  is  oriented  properly,  is  so  small  that 
it  merits  further  examination.  We  will  illustrate  this  by  assuming  a 
buildup  of  slag  in  an  appropriate  orientation  to  produce  a  positive 
yaw  moment  due  to  thrust  and  estimate  the  amount  of  mass  required, 
for  typical  PAM  system  parameters.  We  do  not  attempt  to  justify  the 
assumed  mass  motion,  which  remains  somewhat  an  anomaly. 

Referring  to  Fig.  3,  assume  a  linear  buildup  of  mass  m  =  mt  in 
the  aft  end  of  the  motor  at  r\  =  r  in  the  processing  £tj£  coordinate 
system.  The  axial  acceleration  along  £  due  to  thrust  will  produce 
the  positive  yaw  moment 

Mt  =  ( Tmr/M )  (30) 


The  distance  t  from  the  vehicle  center  of  mass  to  the  nozzle  exit 
plane  where  the  jet-damping  force  acts  is,  in  general,  larger  than  the 
yaw  radius  of  gyration  Jfc,  so  that  the  net  effect  of  the  variable  mass 
is  stabilizing.  Thomson  and  Reiter26  analyzed  variable  mass  effects 
on  the  coning  motion  of  a  solid  rocket  propelled  spacecraft  in  which 
the  vehicle  center  of  mass  displaced  axially  as  the  propellant  burned. 
They  found  the  net  effects  of  variable  mass  to  be  stabilizing,  consis¬ 
tent  with  flight  data.  Other  complex  flowfield  analyses  of  variable 
mass  effects  on  the  PAM  series  of  spin-stabilized  vehicles27*2*  also 
showed  the  net  effects  of  variable  mass  to  be  stabilizing,  which  fails 
to  explain  an  observed  coning  instability  exhibited  by  these  vehi¬ 
cles  toward  the  end  of  the  motor  bum.  The  jet  damping  moment, 
as  defined  in  Refs.  27  ar.d  28,  includes  the  effect  of  variable  yaw 
moment  of  inertia.  It  is  noted  in  Ref.  28  that  their  numerical  predic¬ 
tion  of  jet-damping  moment  is  only  40-80%  of  the  commonly  used 
jet-damping  moment  obtained  from  one-dimensional  flow  theory, 
i.e.,  Eq.  (23).  It  would  appear  from  the  foregoing  analysis  that  this 
reduction  in  jet-damping  moment  is  the  destabilizing  contribution 
from  / //.  The  term  (l2/k2  -  1)  compared  with  i2/k2  is  a  measure 
of  the  reduction  in  jet  damping  due  to  the  variable  inertia. 

Equation  (25)  is  quite  similar  to  Eq.  (27)  in  that  ixflx  has  a 
destabilizing  effect  on  the  roll  rate  whereas  the  jet-damping  term 
is,  in  general,  stabilizing.  If  we  take  the  maximum  value  R  of  the 
ejected  propellant  gases  to  be  the  nozzle  exit  radius  and  assume 
the  rotational  velocity  varies  linearly  with  radius,  like  a  rigid  body 
(per  Ref.  25),  then  F  =  1/2  and  the  jet-damping  term  is  smaller  in 
magnitude  than  ix/Ix  for  typical  PAM  configurations.  Even  if  the 
gas  rotational  velocity  due  to  the  rocket  spin  is  constant  with  radius 
and  equal  to  its  maximum  value  at  the  nozzle  exit  radius,  F  =  2/3 
and  the  destabilizing  effect  of  ix/Ix  still  exceeds  the  jet-damping 
such  that  the  spin  rate  will  increase.  All  of  the  PAM  flight  vehicles 
exhibited  a  small  increase  in  spin  rate  during  the  coning  divergence. 

A  possible  explanation  for  the  PAM  instability  is  the  yaw  mo¬ 
ment  produced  by  the  axial  thrust  acceleration  acting  on  a  small 
mass  of  liquified  slag  that  accumulates  in  the  aft  end  of  the  motor 


which,  when  included  with  variable  inertia  and  jet  damping  in  the 
yaw  equation,  gives 


Tmtr 

ppk2M2 


(31) 


with  the  simplifying  assumptions  of  fixed  center  of  mass  and  con¬ 
stant  yaw  radius  of  gyration.  If  we  assume  that  M  is  a  negative 
constant  such  that  the  vehicle  mass  M  decreases  linearly  from  its 
initial  value  at  ignition  to  its  final  value  at  burnout,  we  can  inte¬ 
grate  Eq.  (31)  for  the  coning  growth  caused  by  the  mass  m.  The 
results  arc  shown  in  Fig.  4  for  only  5  lb  (2.27  kg)  of  mass  buildup  i*» 
the  85-s  bum  time,  with  the  vehicle  properties  of  Table  1  and  with 
t2  fk  2  _  i  selected  to  give  40-80%  of  jet  damping  corresponding  to 
the  results  of  Ref.  28.  The  assumed  slag  buildup  is  only  0.1%  of  the 
burned  propellant  mass.  Although  there  appears  to  be  no  satisfac¬ 
tory  physical  explanation  for  the  mass  motion  prescribed  to  obtain 
these  results,  the  example  illustrates  the  surprisingly  small  amount 
of  mass  required  to  cause  the  coning  divergence  and  the  extreme 
sensitivity  of  the  spinning  vehicle  motion  to  a  yaw  moment 


Table  1  Assumed  PAM/spacecraft 

parameters _ 

Initial  weight  Mog  =  7200  lb  (3266  kg) 
Weight  loss.  Mg  =  -52  Ib/s  (-23.6  kg/s) 

Yaw  radius  of  gyration,  k  —  3  ft  (0.914  m) 

Inertia  ratio,  fx  —  0.4 

Roll  rate,  p  —  300  deg/s 

Thrust  T  =  13, 000  lb  (57,827  N) 

Slag  buildup  rate,  mg  =  5  lb/85  s  (0.0267  kg/s) 
Slag  radial  position,  r  =  1 .5  ft  (0.457  m) 

Initial  coning  angle,  6q  —  0.5  deg 
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Fig.  4  Potential  PAM  coning  divergence  from  5  lb  (2.27  kg)  slag  buildup. 


Missile  Equations  of  Motion 
The  dynamics  of  spinning,  nominally  axisymmetric  missiles  are 
commonly  described  in  terms  of  body-fixed  coordinates29  or  aer- 
oballistic  coordinates  that  pitch  and  yaw  with  the  missile  but  do  not 
roll1-4  Because  the  aerodynamic  forces  and  moments  are  functions 
of  the  missile  angles  of  attack  and  sideslip,  which  depend  on  the 
lateral  translations  of  the  missile  relative  to  its  forward  velocity,  a 
rigorous  formulation  of  the  missile  equations  must  include  the  cou¬ 
pled  translational  and  rotational  equations  of  motion,  Eqs.  (1)  and 
(2).  We  can  obtain  an  approximate  description  of  the  missile  motion 
in  terms  of  the  Euler  angle  moment  equations,  Eqs.  (9—1 1 ),  to  show 
the  similarity  of  spacecraft  and  missile  yaw  moment  instabilities. 

Consider  an  axisymmetric  missile  in  quasisteady  coning  motion 
about  the  mean  flight  path,  as  shown  in  Fig.  5,  and  assume  that 
the  missile  is  statically  stable,  i.e.,  there  is  a  lift  force  L  =  L&9 
proportional  to  9  that  acts  aft  of  the  center  of  mass  a  distance  xs  (the 
static  maigin)  such  to  produce  a  restoring  moment  Lxs  that  tends 
to  reduce  6.  The  missile  center  of  mass  translates  around  the  mean 
flight  path  such  that  the  lift  force,  ignoring  gyroscopic  spin  effects, 
balances  the  centrifugal  force  due  to  precession  according  to 

L  —  m\jr2r  (32) 

where 

r  =  l0  (33) 

if  we  consider  the  missile  rotational  motion  in  the  processing  Euler 
angle  coordinate  frame.  Fig.  1,  and  let  the  missile  axis  lie  along  £, 
then  the  restoring  moment  Lxs  is  —  AT,,,  which  is  defined  in  terms 
of  the  natural  pitch  frequency  (o  —  {L^xJI)  2  according  to 


.2 

myr 


The  length  /  from  the  missile  center  of  mass  to  the  center  of  pre¬ 
cession  can  be  obtained  from  Eqs.  (32)  and  (33)  with  the  result,  Eq. 
(37),  for  the  precession  frequencies,  and  is  found  to  be 

i  =  I /mxs  =  k2/xj  (38) 

Note  that  the  coning  angle  6  in  Fig.  5  is  not  the  missile  total  angle 
of  attack.  Because  the  missile  has  a  lateral  velocity  v  =  ijrr,  as 
shown  in  Fig.  6,  looking  forward  along  the  mean  flight  path,  there  is 
a  component  of  angle  of  attack  normal  to  9  of  magnitude  v/u.  inis 
angle  attack  component  is  proportional  to  0  and  is  found  to  be 


which  is  quite  small,  in  general,  relative  to  9.  The  total  angle  of 
attack  ©  is  then 


If  we  substitute  Eq.  (34)  into  Eq.  (10)  and  assume  9  &  0,  we  can 
write  the  pitch  equation 

(o>2  -f  fipyjr  —  ir2)9  —  0  (35) 

which,  for  nonzero  9 ,  gives  the  two  precession  frequencies 


*+.-  =  Pr  ±  Q  (36) 

For  slender  missiles,  p  <gC  1 ,  and  roll  rates  of  the  order  of  or  less  than 
the  pitch  frequency,  the  precession  frequencies  are  approximately 

=  ±«  (37) 


(40) 


Although  the  lateral  velocity  contribution  to  total  angle  of  attack  is 
generally  quite  small,  it  has  an  appreciable  effect  in  damping  the 
coning  motion,  i.e.,  the  normal  force  or  lift  force  proportional  to 
this  angle-of-attack  component  always  opposes  the  coning  motion 
irrespective  of  whether  \jr  is  positive  or  negative  (clockwise  or  coun¬ 
terclockwise).  Note  also,  from  Fig.  5,  that  a  coning  missile  has  a 
component  of  lateral  rate  $9  about  the  yaw  axis  £.  Aerodynamic 
forces  that  oppose  this  lateral  rate  or  tumbling  motion  also  damp 
the  coning  motion  irrespective  of  its  direction  and  add  to  the  lateral 
force  damping. 
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Fig.  6  Processing  missile  viewed  along  flight  path. 

Yaw  Moment  Undamping 

A  rigorous  derivation  of  the  linearized,  small-angle,  constant  roll 
rate  equations  of  motion  in  terms  of  the  total  angle  of  attack  © 
and  including  aerodynamic  damping  and  a  Magnus  moment  (to  be 
described  later)  gives  for  the  yaw  moment  equation30”32 

(2^  -  pp)&  +  (vjr  -  vmp)e  -  0  (41) 

where  f  has  two  modes,  as  defined  in  Eq.  (36)  or  (37),  and  v  and 
vm  are  defined  in  terms  of  aerodynamic  coefficients  by 


v-ci.+ci.  +  cz. 

(42) 

vm  =  pC*u  +  C*rpa 

(43) 

Equation  (41)  is  of  the  form 

©  +  X+._©  =  0 

(44) 

where  X+t_,  with  \jr  defined  by  Eq.  (36),  is 

* — K-*s) 

(45) 

a  ~vpr  -  vmp  (46) 


and  has  the  solution 

©  =  ©o  exp  (A+>-/)  (47) 

The  motion  described  by  Eqs.  (47)  and  (36),  expressed  in  terms 
of  the  complex  angle  of  attack  f  =  ft  +  ia  in  the  aeroballistic 
coordinate  system,  is  the  damped  epicyclic  motion30-32 

|  =  K+  exp(X+-F  i ir+)t  -I-  AT_  exp(X_  +  i  irj)t  (48) 

where  K+  and  K-  are  determined  from  the  initial  conditions.  The 
parameter  v  in  which  includes  the  lateral  force  plus  lateral  rate 
damping  coefficients,  is  normally  positive  and  damps  both  coning 
motions.  The  parameter  cr,  which  includes  the  Magnus  moment 
coefficient  C*  ,  has  a  damping  effect  on  one  coning  motion  and  an 
undamping  effect  on  the  other.  If  the  roll  rate  and  Magnus  coefficient 
are  large  enough  that  cr/Q  exceeds  v  in  magnitude,  the  exponent 
becomes  positive  for  one  of  the  coning  motions  and  the  coning  angle 
can  grow  to  large  values.  The  Magnus  instability  is  a  classical  case 
of  yaw  moment  undamping  as  illustrated  in  Fig.  7.  The  crossflow 
velocity  due  to  angle  of  attack  adds  to  the  boundary-layer  velocity 
due  to  roll  on  one  side  of  the  vehicle  and  subtracts  from  this  velocity 
on  the  other  side  to  produce  a  net  Magnus  force  or  side  force  normal 
to  the  0  plane.  If  this  force  is  produced  aft  of  the  center  of  mass,  for 
example,  then  the  moment  produced  about  the  center  of  mass  is  the 
positive  yaw  moment  A/f  that  undamps  the  positive  coning  motion 
and  damps  the  negative  coning  motion. 

A  striking  example  of  a  Magnus-type  instability  occurred  with  a 
finned  re-entry  vehicle.33  The  fins  were  canted  to  produce  a  positive 


o 


CROSS  FLOW 
Fig.  7  Magnus  instability. 

(clockwise)  roll  torque.  Because  the  windward  fin  is  more  effective 
than  the  leeward  fin,  with  the  vehicle  at  angle  of  attack,  a  net  side 
force,  or  yaw  force,  is  produced  that  is  approximately  proportional 
to  the  angle  of  attack.  Although  the  mechanism  for  side  force  gen¬ 
eration  is  appreciably  different  from  that  of  the  classical  Magnus 
force  and  is  independent  of  roll  rate  to  the  first  order,  the  net  effect 
is  a  yaw  moment  undamping  that  can  grow  catastrophically,  hence, 
the  designation  “Magnus -type”  instability.  There  are  many  other 
sources  of  side  force  and  Magnus-type  instabilities  with  spinning 
missiles,  including  the  effects  of  re-entry  vehicle  ablation  time  lag,34 
asymmetric  boundary-layer  transition,35  and  vortex  shedding  at  an¬ 
gle  of  attack.36*37  A  particularly  illustrative  example  of  yaw-moment 
undamping  is  roll  resonance  lock-in  from  mass  and  configurational 
asymmetries,  which  is  described  in  the  next  section. 

The  missile  dynamic  stability  is  determined  by  the  sign  of  the 
exponents  in  Eqs.  (47)  and  (48).  For  the  missile  to  remain  dynam¬ 
ically  stable,  the  exponent  must  remain  negative  for  both  preces¬ 
sion  modes.  An  extensive  review  of  symmetric  missile  dynamic 
instabilities  is  contained  in  a  paper  by  Murphy,38  who  considers  lin¬ 
ear  and  nonlinear  Magnus  moments,  nonlinear  static  and  damping 
moments,  configurational  asymmetries,  effects  of  moving  internal 
parts,  roll-induced  side  moments,  unequal  pitch  and  yaw  frequen¬ 
cies,  and  other  aerodynamic  phenomena.  In  most  cases,  dynamic 
instability  can  be  attributed  to  a  yaw  moment  or  to  side  forces  such 
as  the  Magnus  force  of  Fig.  7  that  have  a  component  normal  to  the 
6  plane  that  acts  fore  or  aft  of  the  missile  center  of  mass  to  produce 
a  yaw  moment. 

Roll  Resonance  Lock-In 

Another  classical  example  of  a  yaw  moment  instability  is 
roll  resonance  lock-in  caused  from  mass  and  configurational 
asymmetries.5*39-44  Consider  a  nominally  axisymmetric  missile 
with  small  body-fixed  asymmetries  consisting  of  a  radial  offset  y 
of  the  center  of  mass  from  the  aerodynamic  axis  of  symmetry  and 
an  orthogonal  trim  force  asymmetry  r  that  produces  a  zero  roll  rate 
trim  angle  of  attack  ©r,  as  shown  in  Fig.  8a.  A  quasisteady  solution 
of  the  moment  equations  of  motion  for  the  total  angle  of  attack  0 
yields  the  result32 

-f-  =[(l-A*)*  +  (2fX)*]^  (49) 

©r 

*— "(s)  » 

where  k  —  p/ pent  is  the  ratio  of  the  roll  rate  to  the  critical  roll 
rate,  pent  =  w/(l  /x)i,  and  2f  =  (v  -  P-C*u) /a>(l  -  p)L  Note 
from  Eq.  (50)  and  Fig.  8b  that  for  subcritical  roll  rates  (X  <  1)  4>  is 
positive  and  the  positive  roll  torque  produced  by  the  lift  force  acting 
on  the  moment  arm  y  cos  <p  increases  the  roll  rate  toward  resonance, 
k  —  1.  For  any  positive  value  of  <p  the  trim  component  r  sin  <f>  is  a 
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Fig.  8  Orthogonal  mass  and  configurational  asymmetries. 


yaw  force  that  produces  the  yaw  moment  //  =  w2Qr  sirup  and 
undamps  the  angle  of  attach  As  the  roll  rate  approaches  critical, 
A,  -►  1,  <p  -►  90  deg,  the  trim  angle  of  attack  is  amplified  according 
to  Eq.  (49),  which  rapidly  increases  the  lift  force  and  the  roll  torque 
and  drives  the  missile  into  a  “resonance  lock-in”  condition.  Simi¬ 
larly,  for  initially  supercritical  roll  rates,  X  >  1 ,  <p  is  negative,  which 
produces  a  negative  roll  torque  that  drives  the  roll  rate  down  into 
resonance.  Thus,  the  roll  rate  is  driven  into  the  stable  equilibrium 
condition  of  steady  roll  resonance,  X  =  1,  in  which  the  trim  angle  of 
attack  ©r  is  amplified  by  the  resonance  amplification  factor  1  /2f 
from  Eq.  (49).  The  phenomenon  is  termed  roll  resonance  because 
the  resonance  amplification  of  the  trim  angle  of  attack,  expressed 
by  Eqs.  (49)  and  (50),  is  analogous  to  the  resonance  amplification 
of  a  damped,  single-degree-of-freedom  oscillator.  At  steady  reso¬ 
nance  (<p  =  90  deg)  the  body-fixed  trim  angle  of  attack  produces  a 
pure  yaw  moment,  which  effectively  undamps  the  rolling  trim  angle 
of  attack  coning  motion.  Other  examples  of  roll -resonance-induced 
yaw  moment  instabilities  are  described  in  the  literature.45-49 

Concluding  Remarks 

The  foregoing  derivations  of  the  coning  motions  of  spinning, 
axisymmetric  bodies,  with  selected  examples  of  classical  missile 
and  spacecraft  coning  instabilities,  illustrate  the  extreme  sensitiv¬ 
ity  of  missile  and  spacecraft  motion  to  yaw  moment  disturbances. 
The  analyses  were  performed  in  the  classical  Euler  angle  coordinate 
system,  which  elucidates  the  description  of  the  yaw  moment  and  fa¬ 
cilitates  the  analyses.  The  yaw  moment  is  simply  the  net  moment 
acting  about  one  of  the  coordinate  axes  in  the  precessing  Euler  angle 
system.  It  has  been  shown  that  the  general  motion  of  a  spinning  ax¬ 
isymmetric  body  consists  of  a  precession  about  the  resultant  angular 
momentum  vector  plus  a  spin  of  the  body  about  its  axis  of  symmetry 
relative  to  the  precessing  plane  that  defines  the  coning  angle.  The 
rate  and  direction  of  this  spin  depends  on  the  ratio  of  the  roll  to 
transverse  mass  moment  of  inertia,  which  can  have  a  strong  influ¬ 
ence  on  yaw  moments  generated  by  small  missile  and  spacecraft 
components  free  to  move  relative  to  the  spinning  body.  It  has  been 
shown  that  very  small  asymmetries,  consisting  of  either  mass  or 
configurational  asymmetries  or  combinations  thereof,  can  produce 
small  yaw  moments  sufficient  to  destabilize  the  motion,  sometimes 
leading  to  catastrophic  failures.  It  is  hoped  that  the  viewpoint  ex¬ 
pressed  in  this  paper  has  provided  the  reader  a  bit  more  insight  into 
the  mechanics  of  yaw  moment  undamping  and  can,  perhaps,  help 
the  missile  or  spacecraft  designer  to  avoid  the  occurrence  of  such 
instabilities  in  ftiture  designs. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer”  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation's  Technology  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staffs  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  Micro-Electro-Mechanical  Systems  (MEMS),  and  data  storage 
and  display  technologies;  lasers  and  electro-optics,  solid  state  laser  design,  micro-optics, 
optical  communications,  and  fiber  optic  sensors;  atomic  frequency  standards,  applied 
laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam  control, 
LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of  new 
materials:  metals,  alloys,  ceramics,  polymers  and  composites;  development  and  analysis 
of  advanced  materials  processing  and  deposition  techniques;  nondestructive  evaluation, 
component  failure  analysis  and  reliability;  fracture  mechanics  and  stress  corrosion;  analy¬ 
sis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  spacecraft  structural 
mechanics,  space  environment  effects  on  materials,  hardening  and  vulnerability  assess¬ 
ment;  contamination,  thermal  and  structural  control;  lubrication  and  surface  phenomena; 
microengineering  technology  and  microinstrument  development. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing, 
hyperspectral  imagery;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s  atmos¬ 
phere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate  radia¬ 
tions  on  space  systems;  component  testing,  space  instrumentation;  environmental  moni¬ 
toring,  trace  detection;  atmospheric  chemical  reactions,  atmospheric  optics,  light  scatter¬ 
ing,  state-specific  chemical  reactions  and  radiative  signatures  of  missile  plumes,  and 
sensor  out-of-field-of-view  rejection. 


